Critical Velocities for Energy Dissipation from Periodic Motions of Impurity in 

Bose-Einstein Condensates 
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A phenomenon of energy dissipation in Bose-Einstein condensates is studied based on a micro- 
scopic model for the motion of impurity. Critical velocities for onset of energy dissipation are 
obtained for periodic motions, such as a dipole-like oscillation and a circular motion. The first 
example is similar to a series of MIT group experiments settings where the critical velocity was 
observed much below the Landau critical velocity. The appearance of the smaller values for the crit- 
ical velocity is also found in our model, even in the homogeneous condensate in the thermodynamic 
limit. This suggests that the landau criterion be reexamined in the absence of quantized vortices in 
the bulk limit. 

PACS numbers: 03.75.Fi, 03.75.Kk, 67.40.Yv 



I. INTRODUCTION 

In 1941 Landau gave a phenomenological argument 
about the critical velocity in a superfluid below which 
no energy dissipations occur [1]. His argument was en- 
tirely based on a kinematics and the Galilean invariance. 
This critical velocity is now known as the Landau crite- 
rion, which states that if impurities in a superfluid move 
slower than v c — mino;(p)/|p|, then there are no excita- 
tions created in a superfluid. Here u)(p) is the dispersion 
relation of excitations with the momentum p. This then 
explains the phenomenon of the superfluidity. In a real 
superfluid, however, finite energy dissipations take place 
even below the Landau critical velocity. This is due to 
creation of quantized vortices, which was explained later 
by Feynman [2] . A microscopic understanding of the su- 
perfluidity was initiated by Bogoliubov [3]. He partic- 
ularly showed that the phenomenon of superfluidity can 
be explained as a consequence of condensation of massive 
bosons in the ground state. Although his model cannot 
be strictly speaking, applied to the real 4 He due to strong 
interactions between atoms, he succeed to explain the 
Landau criterion from the first principle. Bogoliubov's 
weakly interacting massive boson model has been stud- 
ied extensively in more than half century [6] . 

The celebrated experimental realization of Bose- 
Einstein condensates (BECs) in several alkali vapors have 
revived theoretical works on BECs as well as experimen- 
tal works [4, 5]. Particularly, the mean field approx- 
imation of the Bogoliubov model, the Gross-Pitaevskii 
equation, has been studied extensively [5]. One reason is 
because this model can be adopted as a realistic model. 
Another reason is that the current experimental technics 
are remarkably well developed and under controlled to 
provide variety of opportunities to test these theories. 
Indeed, so far theories have good agreements with exper- 
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imental results in many cases. However, there are still 
experimental facts which are left without satisfactory ex- 
planation. In this paper we would like to give analytical 
argument to one of those experiments. They are a series 
of experiments done by MIT group [7-9] . 

These experiment were intended to examine the Lan- 
dau criterion in BECs. Interestingly, they have found dis- 
agreement with Landau's argument. The essence of these 
experiments is as follow. The sodium condensate was cre- 
ated in anisotropic harmonic traps forming an elongated 
in one axial direction. A Gaussian laser beam was then 
used to stir the condensate by moving it back and forth 
periodically. The critical velocity for creating energy dis- 
sipation was observed less ten times smaller than the Lan- 
dau critical velocity. Since a publication of these exper- 
imental results, there appear many theoretical works to 
fulfill this discrepancy. One of the important question is 
whether this experimental fact is due to the bulk prop- 
erty of BECs or is originated from other factors, such as 
creation of vortices [10], the geometry of the condensate 
[11], and so on [12]. There is no doubt about the fact 
that all these factors give rise to the experimental obser- 
vations of energy dissipations below the Landau critical 
velocity. However, there still remain unclear whether this 
discrepancy may happen in the homogeneous system. 

In this paper we argue that the appearance of the 
smaller critical velocity is also part of the bulk prop- 
erties of the condensates. Hence this effect does occur 
even in the homogeneous system in the thermodynamic 
limit as contrast to the original Landau's argument. To 
demonstrate it we will evaluate the critical velocity for 
the energy dissipation in the homogeneous BEC in the 
thermodynamics limit at zero temperature. To capture 
the physics of MIT experiments, we couple the weakly 
interacting bosons with an external impurity whose tra- 
jectory is given by a periodic motion similar to the actual 
experiment. We show that this model provides qualita- 
tive understandings of the discrepancy observed in labo- 
ratories. We also compare the result with other possible 
periodic motion, a circular motion. The result shows that 
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the same conclusion holds for the case of the circular mo- 
tion. A comparison between two cases suggests that the 
circular motion has more dissipation than the dipole-like 
oscillation and the larger critical velocity. 

The paper continues as follows. In Sec. II we give a 
summary of our model and its solution within the Bo- 
goliubov approximation. A formula for the energy dis- 
sipation due to the motion of impurity in BECs is also 
given. An idealized case of MIT group experimental set- 
tings is studied and the critical velocity is evaluated in 
Sec. III. In Sec. IV energy dissipation from a circular mo- 
tion and its comparison to the result of Sec. Ill are given. 
Sec. V gives a summary and discussion of our results. 



II. THE MODEL AND ITS SOLUTION 

A. The model Hamiltonian and its diagonalization 

The model for motions of classical impurities in the 
homogeneous BEC was proposed and studied in [18-20]. 
We give a summary of result together with the essence 
of the model. The total Hamiltonian is the sum of two 
parts. One is the standard Bogoliubov weakly interacting 
massive bosons term Hb , and the other is the interaction 
term Hi between bosons and a moving impurity which 
takes into account the local interaction between them : 

Hb= Jd 3 x ft( x ,t)(-^)4(x,t) 

+ | Jd 3 xft(x,t)ft(x,t)j>(x,t)i>(x,t), (1) 

where M is the mass of the bosons. The coupling con- 
stant g between bosons is expressed in terms of the s-wave 
scattering length a s ; g = 4ira s h 2 /M. We assume the re- 
pulsive interaction g > and diluteness na 3 <C 1 with n 
the number density of bosons. The interaction term is 

Hi = X J d 3 x p c (x,t)ft(x,t)ip(x,t). (2) 

Here A is the coupling constant between bosons and the 
impurity, and p c (x,t) is the time dependent distribution 
of classical impurity. Particularly, we have in mind an 
idealized point-like disturbance on the condensate guided 
along a given trajectory £(t), i.e. we will adopt the re- 
placement p c (x, t) = S(x — C(t))- 

In order to estimate the effects of the impurity in 
the thermodynamics limit, we expand the field operators 
xj){x,t) in terms of the plane wave basis with periodic 
boundary conditions in a finite size box V = L 3 . The 
limit N and V go to infinity whit a fixed density n = N/V 
will be taken at the end of calculations. We follow Bo- 
goliubov's treatment to simplify the total Hamiltonian 
H = Hb + Hi within a number conserving framework 
of Bogoliubov model [21]. After several steps [20], we 



arrived at an approximated Hamiltonian : 

H B = E Q + (^k + 9n)&l& k + Y^2 ( & W-k + h - c -) 

+ ( 3 ) 

Here Eq = gnN/2 — gn/2 + \n is a constant term, e k = 
h 2 k 2 /2M is the free kinetic energy of bosons, and p~k(t) = 
J d 3 x p c (x,t)e~ lk ' x . In this approximation we neglected 
terms of order of AT -1 . 

Bogoliubov's excitation is created and annihilated by 
the operators 

= a£cosh0 fe + d_fesinh0 fc , (4) 
bk = dfc cosh0 fc + a { _ k sinh# fe , (5) 

with 

9 k = tanh _1 (- ^— ■ ), hw k = y/e k (e k + 2gn). 

tuJk + e k +gn 

(6) 

In terms of Bogoliubov's excitation, the total Hamilto- 
nian is expressed as 

H = K + E' *o> k i>ib h + Yl W)b{ + h.c). (7) 

E, 
(huj k — e k — gn)/2 is the ground 

state energy without impurities (a prime indicates the 
omission of the zero mode from the summation), and 
fk(t) = n\pk (t) \J e k l (Nhw k ). Therefore, the fist order 
impurity effects are equivalent to decoupled forced har- 
monic oscillators, which can be solved analytically. Since 
we are solving the time dependent problem, it is natu- 
ral to work in the Hcisenberg picture. The equations of 
motions for Bogoliubov's excitation creation and annihi- 
lation operators are easily solved as follows. 

bl(t) = Bl(t) + ct>Ut), (8) 
bk(t) = B k (t) + (f>k(t). (9) 

Here (f>k(t) is a c- number function : 

*^'t(S-^ {t ^ < 10 > 

where the integral Ik(t) is defined by 

/*(*)= [ dt' p k (t')e iuJkt ' ■ (11) 

J to 

We have chosen the boundary condition such that the 
system is disturbed by impurity at t = t . 

Now the our system is described by in terms of the 
dressed Bogoliubov excitation creation and annihilation 
operators B\, and Bk- They evolve according to the di- 
agonalized Hamiltonian H' : 

H' = E (t) + J2' ^kB{B k . (12) 
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In the diagonalized Hamiltonian, the new ground state 

energy is defined by E (t) = E' + 5^Re(/jJ(i)0fc(t)). 
The energy spectrum huj^ is that of the gapless excita- 
tions characterized by u>k — kc s for a small k, where 
c s = y/gn/M is the speed of sound. We remark that 
the spectrum oJk and the speed of sound c are the same 
as in the original Bogoliubov model without impurities. 
Hence, the motion of impurities does not affect either u>k 
nor c in our model within the above approximation. 

B. Homogeneous condensate 

The homogeneous BEC is defined by the ground state 
of the Hamiltonian (7) without the impurity term. In 
other word the vacuum state of the annihilation operators 
fefe. We denote it as \bec), and this is to be defined by 

b k \bec) =0 Vfc ^ 0. (13) 

Therefore, we conclude that the dressed Bogoliubov ex- 
citation annihilation operator due to impurity does not 
see the condensate as vacuum : 

B k \bec) =-<j> k {t)\bec). (14) 

Alternatively, one can interpret this result as follows. A 
classical field describe by <frk{t) is induced in the homo- 
geneous condensate due to the motion of impurity. 

C. Energy dissipation due to the motion of 
impurity 

As we have mentioned in Introduction, we are inter- 
ested in the amount of dissipated energy in condensates 
due to the motion of impurity. Physically, energy dissipa- 
tions take place via scattering processes between bosons 
and impurity. In other words the impurity motion cre- 
ates excitations in condensates by scattering processes. 
Excitations in BECs will carry amount of energy equal 
to their excitation spectra. The amount of dissipated en- 
ergy is then measured by heat transferred to condensates. 

We first define the occupation number for the dressed 
Bogoliubov's excitation with respect to the homogeneous 
condensate. This number counts the emitted Bogoli- 
ubov's excitations accompanying with the motion of im- 
purity in BEC, 

hk(t) = (bec\Bl(t)B k (t)\bec) = |0 fe (t)| 2 (15) 

Multiplying by the excitation energy hu>k gives the dissi- 
pated energy £k(t) for a given mode fc, and the total dis- 
sipated energy £(t) is given by summing over all modes. 
In our model we assume that an impurity is driven by 
some external forces, in which back reaction is negligi- 
ble. Therefore, this dissipated energy is equivalent to 
the amount of energy transfered to the condensate from 
the external impurity. In the thermodynamic limit, the 



total energy dissipation due to impurity is evaluated by 
rather simple formula : 

£{t) = £ # £ fc (t) - «*l'*(*)l a - ( 16 ) 



D. Depletion of the condensate 

The depletion of the condensate d(t) due to the quan- 
tum fluctuation is defined by the formula : 

d(t) = Y^( bec \ & l & k\bec)/N. (17) 

Using the Bogoliubov transformation and the solution 
obtained before, 

+ ^E'(^wi a + ^i«w-^wi% as) 

Since |</>fc(i)| 2 has an additional factor 1/N (see eq. (10)), 
the motion of impurity does not contribute to the de- 
pletion of the homogeneous condensate in the thermo- 
dynamic limit. The depletion is then given by d = 
(8 / '3) y/ na^ / tt which is independent of the impurity ef- 
fects. Therefore, in the large N limit, the homogeneous 
condensate is stable against an external impurity. In 
contrast, however, these terms in the depletion of con- 
densates is not negligible in real experiments where the 
finiteness of the number of particles and finite size effects 
play important roles. 



III. MIT GROUP EXPERIMENTS 

We now apply our model to the experimental settings 
done by MIT group. We consider a point-like impurity 
which is oscillating along the z-axis with a period 2ir/Cl 
and the oscillation amplitude d. The trajectory is ex- 
pressed by C(t) = (0, 0, dcos fit). Although this is still an 
idealized model, we will see that this model gives qual- 
itative explanation for the real experiments. To make 
the corresponding to the experiments, we also define the 
velocity of the impurity by v — 2d/nfl. In order to ob- 
tain analytical expressions for the energy transferred to 
the condensate, we further make an assumption that the 
impurity is oscillating for a sufficiently very long time. 
This assumption seems reasonable if we look the real 
experimental values. Under this assumption we let the 
time integral (11) from — oo to +oo. Using the plane 
wave expansion in terms of the Bessel function 3e(x), 
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exp(— ixcos<fi) — X)fc-oo( — *) £ cxp(i£(/))3g(x), we obtain 



J fc (oo) = ^ (-ij'J/^dcosfl) / dte iuJkt+imt (19) 

fc-oo "'- 00 

OO 

= 27r^(-i) £ J^(/crfcos6i) (5(w fe (20) 



Here the wave number vector fe is written in the spherical 
coordinates with an angle from the z-axis. Therefore, 
in the idealized infinite limit, the integral in considera- 
tion will take discrete values labeled by the integer £. In 
this limit an additional care has to be take to evaluate 
the dissipated energy since we cannot take the square of 
the Dirac delta distributions. To this end we first eval- 
uate the finite time interval T, then we take the limit 
T — > oo in the last step. After straightforward but te- 
dious computations, wc find the total dissipated energy 
per unit time is 



lim 

T^oo 

27rnA 2 



h 2 



gtot(r) 
T 

OO „ 



d 3 k 
(2tt) 3 



ek(3e(kd cos 9)) 5{u>k 



(21) 

m). 

(22) 



The above integrals can be evaluated in terms of the gen- 
eral hypergeometric function : 



2 F 3 (a,b;c,d,e;z) = £ 



n=0 



(o)n(b)n ^ 
{c) n (d)n(e)n Tl\ ' 



(23) 



where (a) n = T(a + n)/T(a) with r(x) the Gamma func- 
tion. The final expression for energy dissipation is 



r = r £ 



{ked/2) 



21 



t{k e £) s 



^ (2£+l)(£!) 2 l + 2(fc/0 2 
3 



x 2 F 3 (£+^£+^£+l 



,2£+l:-(k e d) 2 ). (24) 



Here £ = h/(2Mc s ) is the healing length of the conden- 
sate, and r is a constant with dimension of energy per 
unit time ; T = An\ 2 M 2 c s VL/(irh?'). The wave number 
takes the discrete values as a consequence of the Dirac 
delta distribution as 



kg 



1 



V2£ 



y/1 + {m/c s f 



1 



1/2 



(£ = l,2,---). 



(25) 

We plot the dissipated energy per unit time as a func- 
tion of the velocity of impurity v = 2dn/n. The velocity 
is in units of the speed of sound c s , and the energy per 
unit time is scaled by its value for v = c s . Several curves 
correspond to various values of the oscillation amplitude 
d in units of £; d/£ = 1, 10, 100, and 1000. Except for the 
case d/£ = 1, all curves merge to the same curve. This 
universal behavior shows that critical velocities are the 
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FIG. 1: The dissipated energy per unit time for the dipole- 
like oscillation as a function of the velocity v. The velocity 
is plotted in units of the speed of sound c s , and the energy 
per unit time is scaled by its value for v — c s . There are four 
curves which correspond to the values of amplitude; d/£ = 
1, 10, 100, 1000. 



characteristics of the motion of impurity itself and are 
independent of the parameters involved. This statement 
seems true in general in the bulk limit of condensates. 

We estimate the value for critical velocity v c by fitting 
the curve with the expression T oc v(v — v c ) for small 
values of v. This is based on the standard argument that 
the condensates will experience the drag force from the 
impurity, which is proportional to v — v c . The inner 
product of this force and the velocity of impurity then 
gives the rate of dissipated energy in the condensates. 
The estimated value of the critical velocity for the case 
d/% = 100 is found as v c — 0.34 c s . 



IV. ENERGY DISSIPATION FROM A 
CIRCULAR MOTION 

As an another example for periodic motions, we con- 
sider a circular motion on the xy-plane, which was stud- 
ied in details in [20]. The trajectory is specified by two 
parameters, the radius R and the angular velocity Q; 
C(t) = (RcosClt, iJsinQt, 0). Following the same calcu- 
lation methods carried out for the previous section, the 
formula for the dissipated energy per unit time from the 
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circular motion is 

oo ,,,, oc 

r = E r ° i + (fa0 2 ^ J W« (26) 

where T = n\ 2 MVL / {2tt hi 2 R) , and the discrete wave 
number fe^ is given by eq. (25). Note that in the cir- 
cular motion case, the velocity of the impurity is given 
by v — RQ instead. 

We plot the the dissipated energy per unit time in 
FIG. 2 in the same manner as FIG. 1. The critical veloc- 
ity for the circular motion is estimated from the graph as 
v c = 0.43 c s . This value is slightly bigger than the case 
of the dipole oscillation. 




d 

3)t!}j UOIIEdlSSlQ Agj3ug 



FIG. 2: The dissipated energy per unit time for the circular 
motion as a function of the velocity v. The same convention 
for the units of energy per unit time and the velocity is used 
as in FIG. 1. 

For a comparison, we also plot two cases in FIG. 3 
for the value d/£ = 100 and R/£ = 100. The units 
of the dissipated energy per unit time in FIG. 3 is 
7o = nX 2 /(2h^ 3 ). Two curves show the similarity be- 
tween two cases. Although we cannot make one to one 
corresponding for the velocities of two cases, the circular 
motion has more energy dissipation in general. 

V. CONCLUSIONS AND OUTLOOK 

We have shown that finite energy dissipation from the 
periodic motions of impurity can take place even below 
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FIG. 3: The dissipated energy per unit time for two examples 
for the amplitude d/£ = 100 and the rotation radius R/£ = 
100. The energy per unit time is in units of 70 = nA 2 /(2ft^ 3 ). 



the Landau critical velocity in the homogeneous BEC. 
We have found the critical velocities v c for the dipolc- 
like oscillation and the circular motion as 0.34; c s and 
0.43 c s respectively. Our result qualitatively agrees with 
the experimental observations of Ref. [7-9]. The differ- 
ences arise from the fact that there are many other factors 
as pointed out in previous studies [10-12], as well as sev- 
eral simplifications in our model and calculations. It is, 
therefore, necessary to extend to our model to the case 
for BECs with trapping potentials at finite temperature 
to see whether our model yields better agreement or not. 

The result present suggest that the landau criterion be 
reexamined even in the absence of quantized vortices for 
the homogeneous system. A simple physical argument is 
that the original Landau's argument cannot be applied to 
the case where impurities move under the acceleration. 
This is because the Galilean invariance fails if there is 
an acceleration between two coordinates systems. As a 
related issue, the phenomenon of Cherenkov-like radia- 
tion in BECs was studied for an impurity moving with a 
constant velocity, where the critical velocity for a radia- 
tion was found to be exactly same as the Landau critical 
velocity [16, 19, 20, 25]. 

In this paper we also evaluate the dissipated energy 
from the circular motion of impurity. A trajectory of the 
circular motion in BECs is in principle experimentally 
realizable in current experimental technics. The studies 
of impurities under the circular motion in BECs may 
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give us further understanding of the coherent properties 
of condensates and the phenomenon of superfluidities in 
BECs. 
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